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Abstract 

We present a tutorial on the properties of the new ideal circuit element, a memristor. By 
definition, a memristor M relates the charge q and the magnetic flux in a circuit, and complements 
a resistor R, a capacitor C, and an inductor L as an ingredient of ideal electrical circuits. The 
properties of these three elements and their circuits are a part of the standard curricula. The 
existence of the memristor as the fourth ideal circuit element was predicted in 1971 based on 
symmetry arguments, but was clearly experimentally demonstrated just this year. We present 
the properties of a single memristor, memristors in series and parallel, as well as ideal memristor- 
capacitor (MC), memristor-inductor (ML), and memristor-capacitor-inductor (MCL) circuits. We 
find that the memristor has hysteretic current- voltage characteristics. We show that the ideal MC 
(ML) circuit undergoes non-exponential charge (current) decay with two time-scales, and that by 
switching the polarity of the capacitor, an ideal MCL circuit can be tuned from overdamped to 
underdamped. We present simple models which show that these unusual properties are closely 
related to the memristor's internal dynamics. This tutorial complements the pedagogy of ideal 
circuit elements (R, C, and L) and the properties of their circuits. 
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I. INTRODUCTION 



The properties of basic electrical circuits, constructed from three ideal elements, a resistor, 
a capacitor, an inductor, and an ideal voltage source v(t), are a standard staple of physics 
and engineering courses. These circuits show a wide variety of phenomena such as the 
exponential charging and discharging of a resistor-capacitor (RC) circuit with time constant 
t rc — RC, the exponential rise and decay of the current in a resistor-inductor (RL) circuit 
with time constant trl = L/R, the non-dissipative oscillations in an inductor-capacitor 
(LC) circuit with frequency u LC = 1/yLC, as well as resonant oscillations in a resistor- 
capacitor-inductor (RCL) circuit induced by an alternating-current (AC) voltage source 
with frequency u ~ ujlc- 1 The behavior of these ideal circuits is determined by Kirchoff's 
current law that follows from the continuity equation, and Kirchoff's voltage law. We remind 
the Reader that Kirchoff's voltage law follows from Maxwell's second equation only when 
the time-dependence of the magnetic field created by the current in the circuit is ignored, 
§ E ■ dl = where the line integral of the electric field E is taken over any closed loop in 
the circuit.- The study of elementary circuits with ideal elements provides us with a recipe 
to understand real-world circuits where every capacitor has a finite resistance, every battery 
has an internal resistance, and every resistor has an inductive component; we assume that 
the real-world circuits can be modeled using only the three ideal elements and an ideal 
voltage source. 

An ideal capacitor is defined by the single-valued relationship between the charge q(t) 
and the voltage v(t) via dq = Cdv. Similarly, an ideal resistor is defined by a single- 
valued relationship between the current i(t) and the voltage v(t) via dv = Rdi, and an ideal 
inductor is defined by a single-valued relationship between the magnetic flux <fr(t) and the 
current i(t) via dip = Ldi. These three definitions provide three relations between the four 
fundamental constituents of the circuit theory, namely the charge q, current i, voltage v, 
and magnetic flux <ft (See Figured]). The definition of current, i = dq/dt, and the Lenz's 
law, v = +d(p/dt, give two more relations between the four constituents. (We define the 
flux such that the sign in the Lenz's law is positive). These five relations, shown in Fig. (TJ 
raise a natural question: Why is an element relating the charge q(t) and magnetic flux (p(t) 
missing? Based on this symmetry argument, in 1971 Leon Chua postulated that a new 
ideal element defined by the single- valued relationship d<p = Mdq must exist. He called 
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this element memristor M, a short for memory-resistor.- This ground-breaking hypothesis 
meant that the trio of ideal circuit elements (R,C,L) were not sufficient to model a basic 
real- world circuit (that may have a memristive component as well). In 1976, Leon Chua and 
Sung Kang extended the analysis further to memristive systems.- 1 ^ These seminal articles 
studied the properties of a memristor, the fourth ideal circuit element, and showed that 
diverse systems such as thermistors, Josephson junctions, and ionic transport in neurons, 
described by the Ho dgkins- Huxley model, are special cases of memristive systems.— >^ 

Despite the simplicity and the soundness of the symmetry argument that predicts the 
existence of the fourth ideal element, experimental realization of a quasi-ideal memristor - 
defined by the single-valued relationship dip = Mdq - remained elusive.-^ 1 ^ Early this year, 
Strukov and co-workers^ created, using a nano-scale thin-film device, the first realization of 
a memristor. They presented an elegant physical model in which the memristor is equivalent 
to a time-dependent resistor whose value at time t is linearly proportional to the amount of 
charge q that has passed through it before. This equivalence follows from the memristor's 
definition and Lenz's law, d<j) = Mdq v = M(q)i. It also implies that the memristor value 
- memristance - is measured in the same units as the resistance. 

In this tutorial, we present the properties of basic electrical circuits with a memristor. 
For the most part, this theoretical investigation uses Kirchoff's law and Ohm's law. In 
the next section, we discuss the memristor model presented in Ref. [lC| and analytically 



derive its i-v characteristics. Section II III contains theoretical results for ideal MC and ML 



circuits. We use the linear drift model, presented in Ref. [10, to describe the dependence 
of the effective resistance of the memristor (memristance) on the charge that has passed 
through it. This simplification allows us to obtain analytical closed-form results. We show 
the charge (current) decay "time- const ant" in an ideal MC (ML) circuit depends on the 
polarity of the memristor. Sec. IIVI is intended for advanced students. In this section, we 
present models that characterize the dependence of the memristance on the dopant drift 
inside the memristor. We show that the memristive behavior is amplified when we use 
models that are more realistic than the one used in preceding sections. In Sec. |V]we discuss 
an ideal MCL circuit. We show that depending on the polarity of the memristor, the MCL 
circuit can be overdamped or underdamped, and thus allows far more tunability than an 
ideal RCL circuit. Sec. [VI] concludes the tutorial with a brief discussion. 
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II. A SINGLE MEMRISTOR 

We start this section with the elegant model of a memristor presented in Ref. Q. it 
consisted of a thin film (5 nm thick) with one layer of insulating Ti02 and oxygen-poor 
Ti0 2 -x each, sandwiched between platinum contacts. The oxygen vacancies in the second 
layer behave as charge +2 mobile dopants. These dopants create a doped Ti0 2 region, whose 
resistance is significantly lower than the resistance of the undoped region. The boundary 
between the doped and undoped regions, and therefore the effective resistance of the thin 
film, depends on the position of these dopants. It, in turn, is determined by their mobility 
fin (~ 1CT 10 cm 2 /V.s)^ and the electric field across the doped region.— Figure [2] shows a 
schematic of a memristor of size D (D ~ 10 nm) modeled as two resistors in series, the 
doped region with size w and the undoped region with size (D — w). The effective resistance 
of such a device is 

W f w \ 

M{w) = -ft ON + (l - -^j ftoFF (1) 

where TZon (^lkfi}^ - is the resistance of the memristor if it is completely doped, and 7^off 
is its resistance if it is undoped. Although Eq.(TjQ) is valid for arbitrary values of and 
7?-off, experimentally, the resistance of the doped TiC>2 film is significantly smaller than 
the undoped film, 7^off/^on ~ 10 2 ^> 1 and therefore ATI = (TZopp — 7?on) ~ 7^-off- 
In the presence of a voltage v (t) the current in the memristor is determined by Kirchoff's 
voltage law v (t) = M(w)i(t). The memristive behavior of this system is reflected in the 
time-dependence of size of the doped region w(t). In the simplest model - the linear-drift 
model - the boundary between the doped and the undoped regions drifts at a constant speed 
Vd given by 

-dl =VD = (2) 
where we have used the fact that a current i(t) corresponds to a uniform electric field 
Tl>oi$i{i)/D across the doped region. Since the (oxygen vacancy) dopant drift can either 
expand or contract the doped region, we characterize the "polarity" of a memristor by 
rj — ±1, where rj = +1 corresponds to the expansion of the doped region. We note that 
"switching the memristor polarity" means reversing the battery terminals, or the ± plates 
of a capacitor (in an MC circuit) or reversing the direction of the initial current (in an 
ML circuit). Eqns.fll])-© are used to determine the i-v characteristics of a memristor. 
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Integrating Eq.(|2J) gives 

/.x , /^ON fA Dq{t) 
w[t) = w + r] — q(t) =w + r}——- (3) 

where wq is the initial size of the doped region. Thus, the width of the doped region w(t) 
changes linearly with the amount of charge that has passed through it.— Qo = D 2 / ^ D TZ ON 
is the charge that is required to pass through the memristor for the dopant boundary to 
move through distance D (typical parameters^ imply Qq ~ C). It provides the natural 
scale for charge in a memristive circuit. Substituting this result in Eq.(TjQ) gives 

M(g) = n - V ^, (4) 

vo 

where TZo = TZo^{wq/ D) +7?.off(1 — Wo/D) is the effective resistance (memristance) at time 
t = 0. Eq.dlJ) shows explicitly that the memristance M(q) depends purely on the charge q 
that has passed through it. Combined with v(t) = M(q)i(t), Eq.Q implies that the model 
presented here is an ideal memristor. (We recall that v = M(q)i is equivalent to dip = Mdq). 
The prefactor of the g-dependent term is proportional to 1/D 2 and becomes increasingly 
important when D is small. In addition, for a given D, the memristive effects become 
important only when ATI ^> TZq. Now that we have discussed the memristor model from 



Ref. [10|, in the following paragraphs we obtain analytical results for its i-v characteristics. 

For an ideal circuit with a single memristor and a voltage supply, Kirchoff's voltage law 
implies 

j &nq(t) \ dq d ( AR£\ 

r° - )Tt = it [ Koq - 1 = v[ty (5) 

The solution of this equation, subject to the boundary condition q(0) = is 



ATI 



2ATI 1 



(6) 



m = v{t) , 1 = v{t) (7) 

y/1 - 27 1 ATZ(j ) (t)/QoTil M(q(t)Y 

where <p(t) = j^dTv{r) is the magnetic flux associated with the voltage v(t). Eqs.fjSJ)-© 
provide analytical results for i-v characteristics of an ideal memristor circuit. Eq.Q shows 
that the charge is an invertible function of the magnetic flux 3,4 consistent with the defining 
equation dcj) = M(q)dq. Eq.([7j) shows that a memristor does not introduce a phase-shift 
between the current and the voltage, i = if and only if v = 0. Therefore, unlike an ideal 
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capacitor or an inductor, is a purely dissipative element.- For an AC voltage v(t) = v Q sin (cut), 
the magnetic flux is <f>(t) = vq[1 — cos {wt)} ju) . Note that although v(ir/u — t) = v(t), 
4>(7t/uj — t) 7^ <j){t). Therefore, it follows from Eq.f71) that the current i(v) will be a multi- 
valued function of the voltage v. It also follows that since <fi oc 1/ui, the memristive behavior 
is dominant only at low frequencies uj < uj = 2n/t . Here t = D 2 /fi D v is the time 
that the dopants need to travel distance D under a constant voltage vq. to and ujq provide 
the natural time and frequency scales for a memristive circuit (typical parameters^ imply 
to ~ 0.1 ms and u ~ 50 KHz). We emphasize that Eq.Q is based on the linear-drift 
model, Eq.(j2]), and is valid 12 only when the charge flowing through the memristor is less 
than q max (t) = Qo(l — w /D) when r\ = +1 or q max (t) = QqWq/ D when r\ = — 1. It is easy to 
obtain a diversity of i-v characteristics using Eqns.fjHD an d dZD, including those presented in 



Ref. 
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by choosing appropriate functional forms of v(t). Figure [3] shows the theoretical i-v 
curves for v(t) = v sm(ujt) for u = 0.5uj (red solid), u = u (green dashed), and uj = 5u 
(blue dotted). In each case, the high initial resistance TZq leads to the small slope of the 
i-v curves at the beginning. For uj < ujq as the voltage increases, the size of the doped 
region increases and the memristance decreases. Therefore, the slope of the i-v curve on 
the return sweep is large creating a hysteresis loop. The size of this loop varies inversely 
with the frequency uj. At high frequencies, uj = 5ujq, the size of the doped region barely 
changes before the applied voltage begins the return sweep. Hence the memristance remains 
essentially unchanged and the hysteretic behavior is suppressed. The inset in Fig. |3] shows 
the theoretical q-<j) curve for uj = 0.5a;o that follows from Eq.([6]). 

Thus, a single memristor shows a wide variety of i-v characteristics based on the frequency 
of the applied voltage. Since the mobility of the (oxygen vacancy) dopants is low, memristive 
effects are appreciable only when the memristor size is nano-scale. Now, we consider an ideal 
circuit with two memristors in series (Fig. [2]). It follows from Kirchoff's laws that if two 
memristors Mi and M2 have the same polarity, 771 = 772, they add like regular resistors, 
M(q) = (JZ01 + 7?o2) — r](A7li + /S.lZ2)q{t) / Q whereas when they have opposite polarities, 
V1V2 = — 1, the g-dependent component is suppressed, M(g) = (7?.oi + ^02) — vi^TZi — 
ATZ 2 )q(t)/Q(). The fact that memristors with same polarities add in series leads to the 
possibility of a superlattice of memristors with micron dimensions instead of the nanoscale 
dimensions. We emphasize that a single memristor cannot be scaled up without losing the 
memristive effect because the relative change in the size of the doped region decreases with 
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scaling. A superlattice of nano-scale memristors, on the other hand, will show the same 
memristive effect when scaled up. We leave the problem of two memristors in parallel to 
the Reader. 

These non-trivial properties of an ideal memristor circuit raise the following question: 
What are the properties of basic circuits with a memristor and a capacitor or an inductor? 
(A memristor-resistor circuit is trivial.) We will explore this question in the subsequent 
sections. 



III. IDEAL MC AND ML CIRCUITS 



Let us consider an ideal MC circuit with a capacitor having an initial charge go and 
no voltage source. The effective resistance of the memristor is determined by its polarity 
(whether the doped region increases or decreases), and since the charge decay time-constant 
of the MC circuit depends on its effective resistance, the capacitor discharge will depend on 
the memristor polarity. Kirchoff's voltage law applied to an ideal MC circuit gives 



(8) 



where q(t) is the charge on the capacitor. We emphasize that the g-dependence of the 
memristance here is M c (q) = TZ — f]ATZ(q — q)/Qo because if q is the remaining charge on 
the capacitor, then the charge that has passed through the memristor is (g — <?)• Eq.© 
is integrated by rewriting it as dq/dt = —q/(a + bq) where a = C{1Zq — rjAlZq^/Qo) and 
b = rjCATZ/Qo. We obtain the following implicit equation 



q(t) exp 



r]ATZq{t) 
TZfQo 



g exp 



exp 



(9) 



where Tip = 7Z — r/A7?.go/Qo is the memristance when the entire charge g has passed 
through the memristor .J^ A small t-expansion of Eq.([9]) shows that the initial cur- 
rent i(0) = qo/TZoC is independent of the memristor polarity 77, and the large-t ex- 
pansion shows that the charge on the capacitor decays exponentially, q(t — > 00) = 
g exp(— t/lZpC) exp(r]A7lqo/7lFQo). In the intermediate region, the naive expectation 
Q^) = go exp[— t/M(w(t))C] is not the self-consistent solution of Eq.Q. Therefore, al- 
though a memristor can be thought of as an effective resistor, its effect in an MC circuit is 
not captured by merely substituting its time-dependent value in place of the resistance in 
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an ideal RC circuit. Qualitatively, since the memristance decreases or increases depending 
on its polarity, we expect that when 77 = +1 the MC circuit will discharge faster than an 
RC circuit with same resistance TZq. That RC circuit, in turn, will discharge faster than 
the same MC circuit when r\ = —1. Figure H] shows the theoretical q-t curves obtained 
by (numerically) integrating Eq.(l8"|). These results indeed fulfill our expectations. We note 
that Eq.([9]), obtained using the linear-drift model, is valid for go < Qo(l — wq/D) when 
1] = +1 which guarantees that the final memristance TZf > is always positive.— The 
inset in Fig. H] shows the time-evolution of the size of the doped region w(t) obtained using 
Eq.Q and confirms the applicability of the linear-drift model. We remind the Reader that 
changing the polarity of the memristor can be accomplished by exchanging the ± plates of 
the fully charged capacitor. 

It is now straightforward to understand an ideal MC circuit with a direct-current (DC) 
voltage source Vq and an uncharged capacitor. This problem is the time-reversed version of 
an MC circuit with the capacitor charge go — vqC and no voltage source. The only salient 
difference is that in the present case, the charge passing through the memristor is the same 
as the charge on the capacitor. Using Kirchoff 's voltage law we obtain the following implicit 
result, 



where IZp = 7Z — riA7l(v C)/Q is the memristance when t — > 00. As before, Eq. lfrol) shows 
that when eta = +l(t] = —1), the ideal MC circuit charges faster (slower) than an ideal RC 
circuit with the same resistance TZq. In particular, the capacitor charging time for rj = +1 
(the doped region widens and the memristance reduces with time) decreases steeply as the 
DC voltage vq — ► Qo(l ~wq/D)/C, the maximum voltage at which the linear-drift model is 
applicable.— 

Now we turn our attention to an ML circuit. Ideal RC and RL circuits are described 
by the same differential equation [dqjdt + q/r RC = 0;di/dt + i/trl = 0) with same 
boundary conditions. Therefore they have identical solutions^ q(t) = go exp(— t/r^c*) and 
i(t) = iq exp(— t/TRL). As we will see below, this equivalence breaks down for MC and ML 
circuits. Let us consider an ideal ML circuit with initial current Kirchoff's voltage law 
implies that 




(10) 




(11) 
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The solution of this equation above is given by i(q) = Aq 2 (t) — Bq(t) + i = (q — q + )(q — q- 



where A = r]ATZ/2Q L, B = TZ Q /L, and q± = (Q n /AK) 1 ± y/ 1 - 2r/A^Li /Qo^o 
are the two real roots of i(q) = 0. We integrate the implicit result using partial fractions 
and get 

(12) 



= "ATT 



q + e l / TML — q- 



where t ml = L/1Zq\J\ — 2r]A7lLi /Q 7l 2 ) is characteristic time associated with the ML 
circuit. The current i(t) in the circuit is 

*W = *o 7- , /T ,„ 1 v* ■ (13) 



Eqs.( ll2p - (|13p provide the set of analytical results for an ideal ML circuit. At small-t Eq. ffTS]) 
becomes i(t) = io(l — tlZo/L), whereas the large-t expansion shows that the current decays 
exponentially, i(t — > oo) = i (2Q L/q + A7lT ML ) 2 exp(— t/r ML ). Since t ml depends on the 
polarity of the memristor, tml{v = +1) > t~ml(v = ~ 1); the ML circuit with r\ — +1 dis- 
charges slower than its RL counterpart whereas the same ML circuit with rj = — 1 discharges 
faster than the RL counterpart. Figure [5] shows the theoretical i-t curves for an ML cir- 
cuit obtained from Eq. flTBT) ; these results are consistent with our qualitative analysis. Note 
that the net charge passing through the memristor in an ML circuit is q(t — > oo) = g_(zo). 
Therefore an upper limit on initial current iq for the validity of the linear-drift model^ is 
given by q-(io) < QqWq/D (77 = —1). As in the case of an ideal MC circuit charge, the ML 
circuit current decays steeply as z'o approaches this upper limit. 

Figs. H] and [5] suggest that ideal MC and ML circuits have a one-to-one correspondence 
analogous to the ideal RC and RL circuits. Therefore, it is tempting to think that solution of 
an ideal ML circuit with a DC voltage vq is straightforward. (In a corresponding RL circuit, 
the current asymptotically approaches Vq/R for t ^> trl = L/R). The relevant differential 
equation obtained using Kirchoff's voltage law, 

4 + U- V ^1) ;(() = „„, (M) 



dt 



shows that it is not the case. In an ML circuit, as the current i(t) asymptotically approaches 
its maximum value, it can pump an arbitrarily large charge q(t) = J^i(r)dr through the 
memristor. Hence, for any non-zero voltage, no matter how small, the linear-drift model 
breaks down at large times when w(t) = wq + rjDq(t) / Qq exceeds D (rj — +1) or becomes 
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negative (77 = —1). This failure of the linear-drift model reflects the fact that when the 
(oxygen vacancy) dopants approach either end of the memristor, their drift is strongly 
suppressed by a non-uniform electric field. Thus, unlike the ideal RL circuit, the steady- 
state current in an ideal ML circuit is not solely determined by the resistance TZq but also 
by the inductor. In the following section, we present more realistic models of the dopant 
drift that take into account its suppression near the memristor boundaries. 

IV. MODELS OF NON-LINEAR DOPANT DRIFT 

The linear-drift model used in preceding sections captures the majority of salient features 
of a memristor. It makes the ideal memristor, MC, and ML circuits analytically tractable and 
leads to closed- form results such as Eqs. (j7]), Qj, and f|T3|) . We leave it as an exercise for the 
Reader to verify that these results reduce to their well-known R, RC, and RL counterparts in 
the limit when the memristive effects are negligible, A1Z — * 0. The linear drift model suffers 
from one serious drawback: it does not take into account the boundary effects. Qualitatively, 
the boundary between the doped and undoped regions moves with speed vd in the bulk of 
the memristor, but that speed is strongly suppressed when it approaches either edge, w ~ 
or w ~ D. We modify Eq.Q to reflect this suppression as follows^ 



The window function F(x) satisfies F(0) = F(l) = to ensure no drift at the boundaries. 
The function F(x) is symmetric about x = 1/2 and monotonically increasing over the 
interval < x < 1/2, < F(x) < 1 = F(x = 1/2). These properties guarantee that the 
difference between this model and the linear-drift model, Eq.([2]), vanishes in the bulk of the 
memristor as w — ► D/2. Motivated by this physical picture, we consider a family of window 
functions parameterized by a positive integer p, F p (x) = 1 — (2x — l) 2p . Note that F p (x) 
satisfies all the constraints for any p. The equation F p (x) = has 2 real roots at x = ±1, and 
2(p — 1) complex roots that occur in conjugate pairs. As p increases F p (x) is approximately 
constant over an increasing interval around x = 1/2 and as p — > 00, F p (x) = 1 for all x 
except at x = 0, 1. (For example, 1 — -F p=16 (x) > 0.1 only for x < 0.035 and 1 — x < 0.035.) 
Thus, F p (x) with large p provides an excellent non-linear generalization of the linear-drift 
model without suffering from its limitations. We note that at finite p Eq. (1151) describes a 




(15) 
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memristive systern^^ that is equivalent to an ideal memristor-^ when p — > oo or when the 
linear-drift approximation is applicable. It is instructive to compare the results for large p 
with those for p — 1, F p= i(x) — 4x(l — x), when the window function imposes a non-linear 
drift over the entire region < w < D. 10 For p = 1 it is possible to integrate Eq. (|15p 
analytically and we obtain 

/ \ Dexp (4r]q(t)/Q ) 

w p=1 {q)=w — f — TT. 16 

D + w [exp {4rjq{t)/Q ) - 1\ 

As expected, when the suppression at the boundaries is taken into account, the size of the 
doped region satisfies < w(t) < D for all t and w(t) asymptotically approaches -D(O) when 
rj = +1(— 1). For p > 1, we numerically solve Eq. ffT5l) with Kirchoff's voltage law applied to 
an ideal MCL circuit 

L^ + M(q(t))i(t) + ^- = v(t), (17) 
using the following simple algorithm 

w j+ i = Wj + v — — F [ — )tj,e t (18) 



L C 

Qj+i = qj + ij+itt- (20) 



»r!-%+i)i; > (19) 



Here, e 4 is the discrete time-step and w^, ij and stand for the doped-region width, current, 
and charge at time tj = je t respectively. The algorithm is stable and accurate for small 

t t < lo-V 

Figure [6] compares the theoretical i-v results for a single memristor with two models for 
the dopant drift: a p = 1 model with non- uniform drift over the entire memristor (red solid) 
and a p = 10 model in which the dopant drift is heavily suppressed only near the boundaries 
(green dashed). We see that as p increases, beyond a critical voltage the memristance drops 
rapidly to IZon as the entire memristor is doped. Figure [7] shows theoretical results for a 
discharging ideal MC circuit obtained using two models: one with p — 1 (green dashed for 
7] = +1 and blue dash-dotted for rj = —1) and the other with p = 10 (red solid for rj = +1 
and magenta dotted for rj = —1). The corresponding window functions F p (x) are shown 
in the inset. We observe that the memristive behavior is enhanced as p increases, leading 
to a dramatic difference between the decay times of a single MC circuit when 7] = +1 (red 
solid) and rj = — 1 (magenta dotted). Fig. [7] also shows that fitting the experimental data to 
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these theoretical results can determine the window function that best captures the realistic 
dopant drift for a given sample. 

The properties of ideal MC and ML circuits with an arbitrary voltage are obtained 
by integrating Eqs. fll5|) and f[T7|) using the algorithm described above. However, as the 
discussion in Sec. [I] shows, these circuits significantly differ from their ideal RC and RL 
counterparts only at low frequencies. 

V. OSCILLATIONS AND DAMPING IN AN MCL CIRCUIT 

In this section, we discuss the last remaining elementary circuit, namely an ideal MCL 
circuit. First let us recall the results for an ideal RCL circuit.- For a circuit with no voltage 
source and an initial charge q^, the time-dependent charge on the capacitor is given by 



where u) 2 = oj\ c — (2trl)~ 2 > defines an underdamped circuit and Cj 2 < defines an 
overdamped circuit. The two results are continuous at u = (critically damped circuit). 
Thus, an RCL circuit is tuned through the critical damping when the resistance in the circuit 



The non-linear differential equation describing an MCL circuit, Eq. ljTTl) . cannot be solved 
analytically due to the (/-dependent memristance. Figure [8] shows theoretical q-t curves for 
a single MCL circuit obtained by numerically integrating Eqs. (fT5|) and (|T7|) using p = 50 
window function. When rj = +1 (red solid) the circuit is underdamped because as the 
capacitor discharges the memristance reduces from its initial value TZ . When 77 = — 1 
(dashed green), the discharging capacitor increases the memristance. Therefore, when 77 = 
— 1 the MCL circuit is overdamped. For comparison the blue dotted line shows the theoretical 
q-t result for an ideal RCL circuit with resistance TZq that is chosen such that the circuit is 
close to critically damped, TZq ~ 2\J L/C. Fig. [8] implies that if we exchange the ± plates of 
the capacitor in an MCL circuit, the charge will decays rapidly or oscillate. This property 
is unique to an MCL circuit and arises essentially due to the memristive effects. 

For completeness, we briefly discuss the behavior of an MCL circuit driven by an AC 
voltage source v(t) = v sm(ut), with zero initial charge on the capacitor. For an ideal RCL 
circuit, the steady-state charge q{t) oscillates with the driving frequency u and amplitude 




(21) 



is increased beyond R c = 2^ L/C. 
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Vo/Ly (u 2 — uj 2 lc ) + (lu/trl) 2 . For a given circuit, the maximum amplitude VqVLC / R oc- 
curs at resonance, u = ujlc an d diverges as R — > 0.- Fig. shows theoretical q-t curves for 
an ideal MCL circuit with rj = +1 driven with v(t) = t>o sin(u;o£). The red solid line corre- 
sponds to low LC frequency ulc = 0.1u , the dashed green line corresponds to resonance, 
luq = ujlci an d the dotted blue line corresponds to high LC frequency ujlc — V^^o- We find 
that irrespective of the memristor polarity, the memristive effects are manifest only in the 
transient region. We leave it as an exercise for the Reader to explore the strong transient 
response for u LC < oj and compare it with the steady state response at resonance ujlc = oj. 



VI. DISCUSSION 



In this tutorial, we have presented theoretical properties of the fourth ideal circuit ele- 
ment, the memristor, and of basic circuits that include a memristor. In keeping with the 
revered tradition in physics, the existence of an ideal memristor was predicted in 1971 3 
based purely on symmetry argument^; however, its experimental discovery&i&Eiifi and the 
accompanying elegant physical picture^ 1 ^ took another 37 years. The circuits we discussed 
complement the standard RC, RL, LC, RCL circuits, thus covering all possible circuits that 
can be formed using the four ideal elements (a memristor, a resistor, a capacitor, and an 
inductor) and a voltage source. We have shown in this tutorial that many phenomena - 
the change in the discharge rate of a capacitor when the ± plates are switched or change 
in the current in a circuit when the battery terminals are swapped - are attributable to 
a memristive component in the circuit.— 1 ^ In such real-world circuit can only be 

mapped on to one of the ideal circuits with memristors. 

The primary property of the memristor is the memory of the charge that has passed 
through it, reflected in its effective resistance M(q). Although the microscopic mechanisms 
for this memory can be different,— 1 ^ dimensional analysis implies that the memristor size 
D and mobility provide a unit of magnetic flux D 2 //j,d that characterizes the memristor. 
Although the underlying idea behind a memristor is straightforward, its nano-scale size 
remains the main challenge in creating and experimentally investigating basic electrical 
circuits discussed in this article. 

We conclude this tutorial by mentioning an alternate possibility. It is well-known that an 
RCL circuit is equivalent 1 to a 1-dimensional mass+spring system in which the position y(t) 
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of the point mass is equivalent to the charge q(t), the mass is L, the spring constant is 1/C, 
and the viscous drag force is given by F(v ) = —jv where 7 = R. Therefore, a memristor is 
equivalent to a viscous force with a y-dependent drag coefficient, F M = —j(y)v. Choosing 
l{y) = 7o — ^m/A, where A is the typical stretch of the spring, will create the equivalent 
of an MCL circuit. Since a viscous force naturally occurs in fluids, a vertical mass+spring 
system in which the mass moves inside a fluid with a large vertical viscosity gradient can 
provide a macroscopic realization of the MCL circuit. 
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Figures 




FIG. 1: Relations between four variables of basic electrical circuit theory: the charge q, current 
i, voltage v and the magnetic flux <p. Three well-known ideal circuit elements R, C and L are 
associated with pairs (dv,di), (dq,dv) and (d(j),di) respectively. The top (bottom) horizontal pair 
is related by Lenz's law (definition). This leaves the pair (d<j),dq) unrelated. Leon Chua postulated 
that, due to symmetry, a fourth ideal element (memristor) that relates this pair, d<\> = Mdq, must 
exist. 
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FIG. 2: a: Schematic of a memristor of length D as two resistors in series. The doped region 



(TiO-2-x in Ref. 



1Q ) has resistance 



1Q ) has resistance 1Zq,nw/D and the undoped region (Ti02 in Ref. 
7^-off(1 — w/D). The size of the doped region, with its charge +2 ionic dopants, changes in response 
to the applied voltage and thus alters the effective resistance of the memristor. b: Two memristors 
with the same polarity in series, d and ud represent the doped and undoped regions respectively. In 
this case, the memristive effect is retained because doped regions in both memristors simultaneously 
shrink or expand, c: Two memristors with opposite polarities in series. The net memristive effect 
is suppressed. 
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FIG. 3: Theoretical i-v characteristics of a memristor with applied voltage v(t) = vq sin(u;t) for 
uj = 0.5ujo (red solid), u = ujq (green dashed), and uj = 5o>o (blue dotted). The memristor 
parameters are wq/D = 0.5 and 7^off/^on = 20. The unit of resistance is T^on, the unit of 
voltage is vq, and the unit of current is Iq = Qo/to- We see that the hysteresis is pronounced for 
uj < ujq and suppressed when uj S> ujq. The inset is a typical q-cfi graph showing that the charge q 
is an invertible function of the flux The unit of flux </>o = voto = D 2 /^ld is determined by the 
memristor properties alone (typical parameters^ imply = 10 -2 Wb). 



18 




FIG. 4: Theoretical q-t characteristics of an ideal MC circuit. The memristor parameters are 
wq/D = 0.5 and T^off/T^on = 20. The initial charge on the capacitor is qo/Qo = 0.45 < 
(1 — wq/D) to ensure the validity of linear-drift model^ and C/Cq = 1. The unit of capacitance 
is Co = Qo/vo = to/T^ON- We see that when r\ = +1 (red solid), the capacitor charge in the MC 
circuit decays about twice as fast as when r\ = — 1 (green dashed). The central blue dotted plot 
shows the exponential charge decay of an RC circuit with same initial resistance TZq. The inset 
shows the time-evolution of the boundary between the doped and undoped regions when r/ = +1 
(red solid) and rj = — 1 (green dashed), and confirms that the linear-drift model is valid.— 
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FIG. 5: Theoretical i-t characteristics of an ideal ML circuit. The memristor parameters are 
wq/D = 0.5 and 7^off/7^on = 30. The initial current in the circuit is small, iq/Iq = 0.135, to 
ensure the validity of the linear-drift model 12 that breaks down when io/Io > 0.140, and L/Lq = 30. 
The unit of inductance is Lq = (^q/Iq = to^-ON- We see that when 77 = +1 (red solid), the current 
in the ML circuit decays slower than when rj = — 1 (green dashed) . The central blue dotted plot 
shows the exponential current decay of an RL circuit with same initial resistance IZo. The inset 
shows the time-evolution of the boundary between the doped and undoped regions when 77 = +1 
(red solid) and 77 = —1 (green dashed), and confirms that the linear-drift model is valid.— 
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FIG. 6: Theoretical i-v curves for a memristor with (realistic) dopant drift modeled by window 
functions F p (x) = 1 — (2x — l) 2p with p = 1 (red solid) and p = 10 (green dashed), in the presence 
of an external voltage v(t) = 2vo sm(ujQt/2). The memristor parameters are wq/D = 0.5 and 
T^off/T^on = 50. We see that the memristive behavior is enhanced at p = 10. The slope of the i-v 
curves at small times is the same, TZq , in both cases whereas the slope on return sweep depends 
on the window function. For large p, the return-sweep slope is T^-qn = 1 ^ 1 an d it corresponds 
to a fully doped memristor. 
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FIG. 7: Theoretical q-t curves for an ideal MC circuit with non- linear dopant drift modeled by 
window functions F p (x) with p = 1 and p = 10 shown in the inset. The green dashed (77 = +1) and 
the blue dash-dotted (77 = —1) correspond to p = 1 window function. The red solid (77 = +1) and 
the magenta dotted (77 = — 1) correspond to the p = 10 window function. The horizontal line at 
q/Qo = 0.1 is a guide to the eye. The memristor parameters are wq/D = 0.5 and T^off/T^on = 20. 
The initial charge on the capacitor is qo/Qo = 0.7 and C '/Co = 1. We see that the memristive effect 
is enhanced for large p when 17 = +1. Hence, for large p the two decay time-scales associated with 
77 = +1 (red solid) and r\ = — 1 (magenta dotted) can differ by a factor of TZq/TZon ^ 1- Fitting 
the experimental data to these results can determine the nature of dopant drift in actual samples. 
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FIG. 8: Theoretical q-t curves for an ideal discharging MCL circuit modeled using the window 
function for p = 50. The circuit parameters are wq/D = 0.5, T^off/T^on = 20, L/Lq = 1, 
C/Co = 0.04, and qo/Qo = 2. The initial resistance TZo = 10.5 implies that the corresponding 
ideal RCL circuit, with ujlc = 1/vLC ~ TZq/2L, is close to critically damped. When rj = +1 
(red solid) we see that the MCL circuit is underdamped, whereas when r\ = — 1 (green dashed) it 
is overdamped. Result for the RCL circuit with the same initial resistance TZo is shown in blue 
dotted line. Thus, a single MCL circuit can be driven from overdamped to underdamped behavior 
by simply exchanging the ± plates on the capacitor. 
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FIG. 9: Theoretical q-t curves for an ideal MCL circuit driven by an AC voltage v(t) = vo sin(wo^) 
with r/ = +1. The circuit parameters wq/D = 0.5, T^off/T^-on = 10> L/Lq = 50 are fixed. The 
capacitance is C/C$ = 2 (red solid), C/Cq = 0.02 (green dashed), and C/Cq = 0.01 (blue dotted). 
We see that for ulc < ^o, the amplitude of the transient effects is comparable to the maximum 
amplitude that occurs at resonance, and that the memristive effect disappears in the steady-state 
solution. 
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